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Abstract 

In this paper, we study a routing problem on the Gaussian multiple relay channel, in which nodes employ a decode-and- 
forward coding strategy. We are interested in routes for the information flow through the relays that achieve the highest DF rate. 
We first construct an algorithm that provably finds optimal DF routes. As the algorithm runs in factorial time in the worst case, 
we propose a polynomial time heuristic algorithm that finds an optimal route with high probability. We demonstrate that that the 
optimal (and near optimal) DF routes are good in practice by simulating a distributed DF coding scheme using low density parity 
, | ■ check codes with puncturing and incremental redundancy. 
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I. Introduction 



The wireless channel allows intermediate nodes in a network to overhear the transmissions of a source-destination pair and 
, act as relays. We take an information theoretic approach to understand communication in this network scenario. The single- 
relay channel (SRC) was introduced by van der Meulen [1]. To date, the largest achievable region for the general SRC is 
due to Cover and El Gamal [2], who constructed two coding strategies, commonly referred to as decode-and-forward (DF) 
and compress-and-forward (CF). Chong, Motani, & Garg [3] recently introduced a different decoding technique to give a 
potentially larger achievable region for the general SRC. The SRC was extended to the multiple-relay channel (MRC) by Gupta 
and Kumar [4], and Xie and Kumar [5], who presented an achievable rate region based on DF. The capacity of the MRC is 
not known except for special cases, including the degraded MRC [5] (achievable by DF), the phase fading MRC where the 
relays are within a certain distance from the source [6] (achievable by DF), and mesh networks [7] [8] (achievable by CF). 
In DF, "we can imagine that there is an information flow from the source node... to the destination node..." [5], We call this 
' flow a route. In this paper, we construct an algorithm to find the optimal (i.e., rate maximizing) route for DF on the Gaussian 
CsJ , MRC. However, we note that the algorithm runs in factorial time in the worst case. We exploit properties of the algorithm to 
OS ' desi gn a heuristic algorithm that runs in polynomial time and finds optimal routes with high probability. 

Recently, low-density parity-check (LDPC) codes [9][10][11][12] and Turbo codes [13][14] for the SRC have been con- 
structed based on DF. In this paper, we simulate the DF-based coding strategy using LDPC codes [15] [16] with incremental 
redundancy on the MRC and compare the performance of different routes. 
Our contributions in this paper are as follows. 

1) We construct an algorithm that finds optimal (rate maximizing) routes for DF on the Gaussian MRC. 

2) We construct a heuristic algorithm that runs in polynomial time and finds optimal routes with high probability. 

3) We implement DF on the MRC using LDPC codes [15][16] and demonstrate that the DF optimal routes are good in 
' practice. 

II. Network Model 

We consider a D-node Gaussian MRC: S = {1, 2,3..., D — 1, D} with one source (node 1) and one destination (node D). 
We use the stand ard pa th loss model for signal propagation. The received signal at node t, t = 2, 3, . . . , D, can be written as 

Yt = Y^x\jtt \J K d-i t r 'Xi + Zt, where the transmit signal of node i, X%, is a random variable with average power constraint 
E[Xf] < d it is the distance between nodes i and t, r\ is the path loss exponent (77 > 2 with equality for free space 
transmission), and k is a positive constant. Z t is the receiver noise at node t, which is an independent, zero mean Gaussian 
random variable with variance iV t . We assume duplex nodes and that all nodes have the same noise variance. 

A. Route 

Now, we define what we mean by a route in a network. 

Definition 1: The route taken by a packet from the source to the destination is an ordered set of nodes involved in 
encoding/transmitting of the packet. The sequence of the nodes in the route is determined by the order in which nodes' 
transmit signals first depend on the packet. Note that the destination node is the last node in the route and does not transmit 
the data. 

Remark 1: For DF, a route is also the order in which a message is fully decoded at the nodes. 



In the rest of this paper, we denote a route by Ai — {mi, m-2, ■ ■ ■ , mi_M\}. We define the set of all possible routes from the 
source (node 1) to the destination (node D) by 11(5) = | {mi, 7712, . . . ,m\M\} '■ m 2j ■ ■ ■ , m \M\-i are all possible selections 

and permutations of the relays (including the empty set), mi = 1, m\M\ = Dj, 

It has been noted [5] [6] that the DF rate depends on the route selected (known as node order/permutation in the papers). 
Besides the degraded case, finding the optimal route is not easy in general. We refer to the strategy of testing all possible 
routes as brute force. 

B. Achievable Rates using DF 

Using DF [5][6][17] with Gaussian inputs and route Ai, node rrii transmits X mi = Y^pll+i \f a nnmj Prm U mj , for < 
Ysj^i+i a mimj < 1, Vi = 1, . . . , \Ai\ — 1. U mj are independent Gaussian random variables with unit variance. {a mirrij \ j = 
. . . , \Ai\} are the power splits of node rrii, allocating portions of its transmit power to transmit independent sub-codewords 

U mj . 

On route Ai, DF can achieve rates up to 

Ruf(M) — max min R mt (M), (1) 

OLij rn t eM\{mi} 

where R mt (Ai) is the reception rate of the node m t in route Ai, given by 




R mt (M) = L \N- > I > j Ja mimj P mimt , (2) 



where P mjTOt = «rf m ^ m P mj and L{x) — |log(l + x). We call the rate in (Q} the DF rate supported by route M.. The 
maximum DF rate is 

i?™ x = max Rap(M). (3) 

C. The Optimal Routing Problem 

We define the optimal route set for DF as 

Qdf = {M £ U(S)\R DF (M) = R^ x }, 

We define the optimal route set because the rate maximizing route may not be unique. Then the optimal DF routing problem 
for a network S is to find A4^p £ Qdf- 

III. The Nearest Neighbor Algorithm 

In this section, we present an algorithm to find an optimal DF route. First, we define the nearest neighbor with respect to 
a route. 

Definition 2: Node i ^ Ai is a nearest neighbor with respect to the route Ai iff 

Pmi > Pmj, Vm e M, Vj G S \ (M U {»}). (4) 
Now, we describe the nearest neighbor algorithm (NNA). 
Algorithm 1 (NNA): 

1) First, start with the source node, Ai — {mi}. 

2) If there exists a unique nearest neighbor i* with respect to the current route Ai, we append i* to the current route: 
A'l ^ U{i*}. Else, the NNA terminates prematurely. Since Ai is an ordered set, the notation AUB means appending 
ordered set B to the end of ordered set A. 

3) Step |2] is repeated until the destination node, node D, is added into Ai. 

The algorithm is said to terminate normally if node D is added to the route. Otherwise, the algorithm is said to terminate 
prematurely. If the NNA terminates normally, we have the following theorem. 

Theorem 1: Consider a Gaussian MRC. If the NNA terminates normally, then the NNA route is optimal for DF. 
To prove Theorem Q] we need the following lemmas. 

Lemma 1: When we add the unique nearest neighbor, node a*, to route Ai, the rate supported by the new route Ai\ — 
Ai U {mixi+1 = a *} is greater than or equal to the rate supported by the route formed by adding any other node to Ai, 
M2 = M U {m\M\+i = b}- Mathematically, 

R DF {Ai U {a*}) > R DF {M U {b}), V6 e S \ (Ai U {a*}). (5) 



Proof: [Proof for Lemma [T] Considering A4 2 , the reception rate at node 





b is 



Here, we choose the set of {a\j } that maximizes min 2 
of the node m\M\+i — a * m route M\ is 

Ra*(Mi) = L (n-* 1 

(2) 

Note that this set of {a\j } does not necessarily maximizes min 2 <i<| J vi|+i RmA-M-i)- Clearly, if P ma * > P m b, Vm G M. 
with at least one inequality and N a * = Nb, we get R a *(Mi) > R b (M.2)- 

Using {o^}, the reception rates of the first \M\ nodes remain the same for both A4i and M. 2 . Using the power split 
optimized for M.\, we might be able to increase the minimum reception rate in M.\. Hence R DF (A4i) > Rdf(M 2 ). ■ 

Next, we show that choosing the nearest neighbor will not harm the rate supported by the route even when more nodes are 
added. 

Lemma 2: Let M. = {a*, a|, • ■ • , a*^} be a route formed by adding the nearest neighbor (assuming it exists) one by one 
starting from the source. Now, arbitrarily add K nodes to M.. The first node b\ is not a nearest neighbor and the rest may or 



(6) 



<i<\M\+i RmA-M-2)- Using the same set of {a\ 3 }, the reception rate 




(7) 



may not be nearest neighbors. In other words, Ai\ = {a{, aj, 



l \M\ 



, b\, 62, . . • , bx }, where b\ is not a nearest neighbor to 



M. We can always choose the nearest neighbor a* M \ +1 (assuming it exists) at the (\M\ + l)-th position, i.e., 

{a*, . . . , a* M \,a* M \ +1 ,bi, . . . , b K -i}, 



M 2 = 



if ajWi ^ {bi,...,b K -i}, 



|M|' a \M\ 



b\, . . . , bk~i,bk+i, 
if a |A4|+i = b k, for some b k e i b i, ■ 



,b K -i} 



(8) 



and show that i?DF(A4 2 ) > Rdf(Mi). 

Proof: [Proof for Lemma[2) First, we study the case when a* M ^ +1 ^ {bi, 

for A4i be {a^}. Using this set {ck^} for the nodes in M 2 and some {a' a 



, bx-i}- Let the optimum set of power splits 
■ Mi+ij }fox node aJ M]+v 

R a j Ml+1 (M 2 )>R bl (M 1 ) (9a) 
Ra?(M 2 ) = R a *(Mi), Vi = 2,3,...,|A4| (9b) 
Rb l {M 2 )>R bl {M 1 ), Vi = l,2,...,K-l. (9c) 

Equation d9at is due to Lemma [T] Equation ( |9bl is because nodes transmitting to a* in both routes uses the same power split. 
Equation ( |9c| > is because all the nodes except a* M \+i uses the same power split for both routes and with an additional node 
a *Ai|+i transmitting to nodes {bi, . . . , &a--i} in M-i, there is a possible increase in the reception rate. 

Now, we study the case when a|^| +1 = bk for some bk G . . . , bx-i}- Now, using {o^} for all the nodes in A4 2 , 

2,3,..., \M\ (10a) 



Ra?(M 2 ) = R a ;(Mi), Vz 
Raj Ml+1 (M 2 ) > R bl (Mi), 

R bi (M 2 ) = R bi {Mi), Vi = l, 



,k- l,fc + l, 



(10b) 
(10c) 



Equations (11 Oat and (11 Oct are because nodes transmit using the same power split {a-, } in both routes. Equation fllObt follows 
from Lemma [T] 

, to 



Now, let {o^ 7 } be a set of power splits optimized for A4 2 . 

Rdf{M 2 ) 

= , min R k (M 2 ), 

keM 2 \{al} 

> min Rk(M 2 ), 
keM 2 \{ai} 

> min R k (Mi), 

keMi\{al} 

= R DF (Mi). 



< ( 2 )l 

using {al/} 



using {a' /}, {a' a , } 

L J " J \M\ + 1 J 



(11a) 
(lib) 

(11c) 

(lid) 

(lie) 



This proves Lemma [2] ■ 
Lemma 3: The DF rate supported by a route that contains all nearest neighbors is always higher or equal to that by any 

route of the same length, with non-nearest neighbor(s). 

Proof: [Proof for Lemma [3] Lemma [3] can be proven by applying Lemma [2] recursively until the entire route is replaced 

by nearest neighbor nodes. ■ 
Now we consider routes from the source to the destination, with possibly different lengths. 

Proof: [Proof for Theorem [T] Consider a route Mi from the source to the destination, which contains one or more 

non-nearest-neighbor(s). If NNA terminates normally and outputs M 2 , we want to show that Rdf(M 2 ) > -Rdf(TWi). We note 

that | Mi | does not necessarily equal |M 2 |. 

Let M\ = {ml = 1, m 2 , . . . , m\ Ml \ = D} and M 2 = {ml — 1, m 2 , . . . , m f M i = D}. We use asterisks to mark nearest 

neighbor nodes. 

First of all, we consider the case |Mi| = |M 2 |. The results follows immediately from Lemma [3] Second, we consider 
|Mi| > |M 2 |. We consider first \M 2 \ nodes in Mi, i.e., M[ — {m\,m 2 , ■ ■ ■ ,W|m,|}- Then, 

Rdf(M 2 ) > Rdf(M[) > Rdf(Mi). (12) 

The first inequality is obtained by applying Lemma|3] as \M 2 \ = \M[\. The second inequality can be argued as follows. The 
first \M 2 \ nodes in both routes M[ and Mi are identical. Hence the reception rates are the same using the same power split. 
However, there are additional nodes in Mi whose reception rate might be lower than Ruf(M[). Hence, Rdf(Mi) < Rdf{M[). 
Lastly, consider |M 2 > Mi|. We replace the transmitting nodes in Mi with nearest neighbors and obtain 

M 3 = {ml,m 2 , . . ^mj^i^miA^i = D}. (13) 

Note that m\M t \ is n °t be the nearest neighbor. Clearly, using Lemma [2] Rf>f(M 3 ) > Rdf(Mi). 
Now, using the set of power splits {o^} optimized for M3, we have 

Rdf{M 2 ) (14a) 
= min R k (M 2 ), using {a\f} (14b) 

k£M 2 \{al} J 

min Rk(M 2 ), using {a,-^} and some {a^} for 

kEM 2 \{al} J 

nodes mJ Mi{ ,. .. ,mJ M2hl (14c) 
> min R k (M 3 ), using {a\f} (14d) 

keM 3 \{al} J 

= Rdf(M 3 ) > Rdf{Mi). (14e) 

The inequality in (114-db is because using the NNA, {mJ Mi ^, . . . , m* M ^ _ 1 } are added to {m*, . . . , b e f° re D. A 

necessary condition for this is 

P mn > P m D,Vm 6 {ml, . . ^mj^i-i^Vn e {m* Mil , . . . , m*^^}, (15) 

Hence, 

R n {M 2 ) > R D (M 3 ),Vn e {mf JMl | J ... > mf jMa ,_ 1 }. (16) 

With additional nodes transmitting to D in M 2 , Rd(M 2 ) > Rd{M 3 ). Hence, we have TheoremQ] ■ 
Remark 2: We note the NNA terminates normally if and only if a unique nearest neighbor exists at each step. In the next 
section, we extend the NNA to an algorithm which terminates normally given any network topology. 

IV. The Nearest Neighbor Set Algorithm 

In this section, we modify the NNA so that it terminates normally in any Gaussian MRC. We term this algorithm the nearest 
neighbor set algorithm (NNSA). First, we define the nearest neighbor set. 

Definition 3: The nearest neighbor set TV = {m, n 2 , . . . , n|jv|} with respect to route M = {mi, m 2 , . . . , is defined 

as the smallest set N where each n 6 J\f C S \ M satisfies the following condition. 

P mn > P ma , Vm e M,VaeS\ (Mil AT), (17) 

with at least one strict inequality for every pair of (n, a) G {(n, a)\n G TV, a G S \ (M U TV)}. 
Now we describe the NNSA. 
Algorithm 2 (NNSA): 

1) Starting with the source node, we have M = {!}. 



2) Find the nearest neighbor set N. The original route M. branches out to \M\ new routes as follows: 



M l ^MU{n l }, i = l,...,\M\. (18) 
3) For each new route in ( fT~8b . step [2] is repeated until the destination is added to all routes. 

When the algorithm terminates, we end up with many routes from the source to the destination. We term these routes NNSA 
candidates. We calculate the supported rate of each candidate and choose the one which gives the highest supported rate. The 
following theorem says that the NNSA candidate that gives the highest supported rate is an optimal route for DF. 

Theorem 2: Consider a Gaussian MRC. The NNSA candidate routes that give the highest supported rate are optimal for 
DF. 

Proof: [Sketch of proof for Theorem |3 Using the technique used in the proof of Theorem Q] we can show that adding a 
node that does not belong to the nearest neighbor set can only be suboptimal. We can always replace that node with one from 
the nearest neighbor set and obtain an equal or higher rate. In other words, consider a non-NNSA candidate Aii = {ml = 
1, m2, ■ ■ ■ i W|Mi| = D}, where one or more nodes in {7712, . . . , 777|Mi|} are not from the nearest neighbor set, and an NNSA 
candidate M.2 = {m* = 1, m^, ■ ■ ■ , m \M 2 \ = D}, where all nodes in M.2 are added according to the NNSA. We can show 
that R DF {M 2 ) > Rdf{Mi). 

The NNSA finds all possible routes for which every node is added from the nearest neighbor set. Hence one or more of the 
NNSA candidate must achieve the highest DF rate. This gives us Theorem [2] ■ 

Remark 3: We can show that a shortest optimal route, A^^p 11 E Qdf, s.t. A^ D p R < |.M|,V./V( G Qdf, is contained in one 
of the NNSA candidates that supports i?g F x . 

Remark 4: We note that the NNSA is also optimum in the phase fading Gaussian MRC, where all node transmit independent 
codewords, i.e., a,j = l,Vi € M. \ {D },j = i+ 1 and ay = 0,Vj + 

V. Complexity of the NNSA 

With the NNSA, we can now search for the optimal route in the NNSA candidate set, as compared to searching in 11(5) 
using brute force. The number of candidates determines the number of routes whose rate we need to calculate to find optimal 
routes. We note that the size of the NNSA candidate set might still, in the worst case, equal |II(<S)|. Using brute force, the 
number of permutations we need to check is 

where D is the total number of nodes in the network and (™) = -, — " x ("~ 1 ) >< '' xl — - 

\kJ in— k) X (n — k — 1) X ■■■ X 1 

We ran the NNSA on 10000 randomly generated networks with a varying number of nodes uniformly distributed in a 
lmx lm square area. The source, relays, and the destination were randomly assigned. On average, half of the NNSA candidate 
set sizes were less than 0.715% of |TT(<S) | for the 8-node channel and less than 0.253% of |TI(<S) | for the 11-node channel. 

We note that the average size of the NNSA candidate set does grow factorially with the number of nodes. However this does 
increase the range of finite size networks for which we can find optimal routes. Furthermore, the NNSA provides insights for 
designing heuristic algorithms to find good routes for DF-based codes. In the next section, we propose a heuristic algorithm 
which outputs routes in polynomial time. 

VI. A Heuristic Algorithm 

In the NNSA, the optimal route is constructed by adding the "next hop" node one by one to the partial route. The node to 
be added is from the nearest neighbor set. If the nearest neighbor set contains more than one node, the current route branches 
to more than one routes, leading to a possibly large NNSA candidate set size. 

To avoid this, we consider a heuristic approach that starts from the source node and repeatedly adds only one "good" 
candidate from the nearest neighbor set until the destination is reached. For the choice of the next hop node, we consider 
the node which receives the largest sum of received power from all the nodes in the existing partial route. We call this the 
maximum sum-of-received-power algorithm (MSPA). By choosing only one node to be added to the partial route, we prevent 
the algorithm from branching out to multiple routes. This heuristic approach yields only one route, regardless of the network 
size. We now explicitly describe the MSPA. 

Algorithm 3 (MSPA): 

1) First, start with the source node, Ai — {mi}. 

2) For every node t € S \ M., find the sum of received power from all nodes in M. to t, 2~2i£M -^»* - 

3) Let a* be any node with the highest sum of received power, i.e., 2~2ieM ^ ia ' — 2~2jeM Pjt^t £5 \ A4. Append node 
a* to the route: M <— M U {a*}. 

4) Repeat steps [2|-[3] until the destination is added to the route. 



TABLE I 
Performance of the MSPA. 



D 


Average R D F (Mmspa) / R'S™ 


Fraction of optimal routes using MSPA 


3 


1 


1 


4 


0.9999950 


0.99882 


5 


0.9999513 


0.99522 


6 


0.9999399 


0.99194 



LDPC in 4-node MRC. 1(0,0.5) 2(0,0.4) 3(0,1) 4(0,1.5) 




(a) A 4-node network (b) Information BER versus transmit SNR 

Fig. 1, Performance of the LDPC coded DF strategy for different routes in the 4-node network. 



Remark 5: Assuming that the value of the previous sum-of-received-power computations are cached, the complexity of 
step|2]in MSPA is O(D) because there are at most (D — 1) nodes not in the route. The complexity of the comparisons in step 
[3]is O(D). Steps |2H3] are repeated at most (D — 1) times, giving a worst case complexity of the MSPA of 0(D 2 ). Recall that 
D = \S\. 

We can show that the MSPA is optimal for DF on the phase fading Gaussian MRC. 

Theorem 3: In a Gaussian MRC in which the nodes send independent codewords, the MSPA route is optimal for DF. 
Proof: [Proof for Theorem [3] See [18] for proof. ■ 

We now investigate how well the MSPA route .Mmspa performs compared to the optimal route in the Gaussian MRC. Due 
to the complexity involved in optimizing the power splits, we only simulate MRC up to 6 nodes, i.e., D < 6. For each D, we 
randomly place D nodes in a (D — l)m x (D — l)m square area. The source and the destination are randomly chosen. We 
run the NNSA to find the optimal rate R^™, and run the MSPA to find -Rdf (.Mmspa)- The results are shown in Table U With 
high probability, the MSPA is able to find an optimal route. Also, i?DF (-Mmspa) is a good indicator of i?pp x - 

VII. DF with LDPC Codes 

In the previous section, we computed achievable rates of different routes using an information theoretic approach. In this 
section, we present the performance (via simulation) in a network using practical LDPC codes [15] [16] with incremental 
redundancy. In the literature, codes have been designed for the SRC [9] [10] [11] [12]. Here, we implement a DF coding strategy 
on the MRC using rate-i LDPC codes with puncturing and incremental redundancy. 

Fig. Q] gives the bit error rate (BER) versus transmit signal-to-noise ratio (SNR) for different routes. We find that different 
routes have significantly different performance (up to 4dB at BER=0.001 in this example). We also find that the NNSA 
route (also the MSPA route), i.e., .Mmspa = -Mdf = {1,2,3,4}, gives the best performance. Furthermore, in this setup, the 
NNSA/MSPA yields the optimal route directly, meaning that we do not need to evaluate all possible routes (5 routes in this 
example). 

Remark 6: In practice, a relay in the route might decode a message wrongly and thus forward the wrong message. When 
this happens, the nodes behind, when trying to cancel the co-channel interference introduced by this relay, will introduce more 
noise at their decoders. While this scenario is not captured in (Q~|), we allow imperfect interference cancellation in all of our 
simulations here. 
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